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Abstract 

A recurrence relation between equal mass two-loop sunrise diagrams differing in dimensionality 
by 2 is derived and it's solution in terms of Gauss' 2F1 and Appell's F2 hypergeometric functions 
is presented. For arbitrary space-time dimension d the imaginary part of the diagram on the cut is 
found to be the 2F1 hypergeometric function with argument proportional to the maximum of the 
Kibble cubic form. The analytic expression for the threshold value of the diagram in terms of the 
hypergeometric function 3_F2 of argument —1/3 is given. 
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1 Introduction 



The evaluation of radiative corrections for modern high precision particle physics is becoming a more 
and more demanding task. Without inventing new mathematical methods and new computer al- 
gorithms the progress in calculating multi-loop, multi-leg Feynman diagrams depending on several 
momentum and mass scales will be not possible. 

An important class of radiative corrections comes from self-energy type of Feynman diagrams, 
which also occur in evaluating vertex, box and higher multileg diagrams. At the two-loop level different 
approaches for calculating self-energy diagrams are available [5- A general algorithm for reduction of 
propagator type of diagrams to a minimal set of master integrals was proposed in |2j. This recurrence 
relations algorithm has been implemented in computer packages in |3] and At present the 

most advanced package available for calculating two-loop self-energy diagrams with arbitrary massive 
particles was written by Stephen Martin and David Robertson [S] . It includes procedures for numerical 
evaluation of master integrals with arbitrary masses and also a database of analytically known master 
integrals. Integral representation for master integrals with arbitrary masses in four dimensional space 
time was proposed in 

Despite intensive efforts by many authors not all two-loop self-energy integrals with a mass are 
known analytically. Even the imaginary part of the simplest sunrise self-energy diagram with three 
equal mass propagators was not known for arbitrary space-time dimension d until now. The two-loop 
sunrise integral with three equal masses was investigated in many publications [Zj- jT^]. Small and 
large momentum expansions of this integral for arbitrary space-time dimension d can be found in 
It's threshold expansion was given in ^M,- A numerical procedure for evaluating the sunrise integral 
was described in The very latest effort of an analytic calculation of the diagram, by using the 
differential equation approach was undertaken in Ref. [f7j. 

It is the purpose of this paper to describe a new method and to present the analytic result for the 
equal mass two-loop sunrise master integral. To accomplish our goal we use the method of evaluation 
of master integrals by dimensional recurrences proposed in [^. The application of this method to 
one-loop integrals was presented in [12] and jlQ| . In the present publication we extend the method to 
two-loop integrals. 

As was already discovered in the one-loop case, the solutions of dimensional recurrences are com- 
binations of hypergeometric functions. The knowledge of the hypergeometric representation of an 
integral means that we possess the most complete mathematical information available. This infor- 
mation can be effectively used in several respects. First, through analytic continuation formulae, the 
hypergeometric functions valid in one kinematic domain can be re-expressed in a different kinematic 
region. Second, these hypergeometric functions often have integral representation themselves, in which 
an expansion in e = (4 — d)/2 can be made, yielding expressions in logarithms, dilogarithms, elliptic 
integrals, etc.. Since very similar hypergeometric functions come from different kind of Feynman in- 
tegrals the e expansion derived in solving one problem can be used in other applications. Essential 
progress in the e expansion of hypergeometric functions encountered in evaluating Feynman diagrams 
was achieved in jllj.jl3j. |21j and [^. Third, because the hypergeometric series is convergent and 
well behaved in a particular region of kinematical variables, it can be numerically evaluated |2S], |24j . 
In addition a hypergeometric representation allows an asymptotic expansion of the integral in terms 
of ratios of different Gram determinants or ratios of momentum and mass scales which can provide 
fast numerical convergence of the result. 

Our paper is organized as follows. In Sec. 2, we present the relevant difference equation connecting 
sunrise integrals with dimensionality differing by 2 as well as the differential equation for this integral. 
In Sec. 3. the method for finding the full solution of the dimensional recurrency is elaborated. Explicit 
expression for the sunrise integral in terms of Appell's function F2 and Gauss' hypergeometric function 
2F1 is constructed in Sec. 4 and in Sec. 5 the differential equation approach and the method of 
dimensional recurrences is compared. In the Appendix some useful formulae for the hypergeometric 
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functions 2-^1 ; -^1 and F2 are given together with their integral representations. 



2 Difference and differential equations for the sunrise integral 



The generic two-loop self-energy type diagram in d dimensional Minkowski space with three equal 
mass propagators is given by the following integral: 

d'^kid'^k2 1 

{md/^Y {kf - m2)^'i {{ki - - m?Y^ {{k2 - qY - m^Y^ ' 
For integer values of Uj the integrals H2.1() can be expressed in terms of only three basis integrals 
4^\l, 1, 1), jf\2, 1, 1) and 1, 1) = (rf )(m2))2 where 

^id),2._ f d'^k 1 _ / d\ ,2 



TY'(m') = / ; , , , = - r 1 - - m^-^ (2.2 



The relation connecting d — 2 and d dimensional integrals J^^\iji, 1^2,1^3) follows from the relationship 
given in Ref. [2]: 

j'^^~'^\vi,V2: 1^3) = l^lV2j'f\vi + 1,1^2 + 1, U3) 

+ UiU3ji^\ui + l,U2,U3 + l)+'^2iy3A'^\l^l,iy2 + l,1^3 + l)- (2.3) 

Relation (|2.H|) taken at j^i = 2^2 = 1^3 = 1 and 1^1 = 2, 1^2 = ^^3 = 1 gives two equations. To simplify 
these equations we use the recurrence relations proposed in j^. From these two equations by shifting 
d — > d + 2 two more relations follow. They are used to exclude j'^\2, 1, 1) from one of the relations, 
so that we obtain a difference equation for the master integral 1, 1) = ^3'^'': 

I2z^(d + l)(d - l)(3ti + 4) (3d + 2) J^^^^'^ 



-4m^(l -3z)(l -42z + 9z^)z(d- l)d 
-Am\l- zY{l-9zY 

3z[{z + l){27z'^ + 18z -l)d'^ -4z{l + 9z)d-A8z^]m^'^+^ , (2.4) 



d-^ 



where 

2 

m 

z = —. 



(2.5) 



The integral J^f^ satisfies also a second order differential equation 11 . Taking the second derivative 



of 4"^^ with respect to mass gives 



1, 1) = 6jf (2, 2, 1) + 6jf (3, 1, 1). (2.6) 



Again using the recurrence relations from [2], the integrals on the r.h.s can be reduced to the same 
three basis integrals. Using 

4'^{2,l,l) = l^4'\l,l,l) (2.7) 



from 1)2. 6|) we obtain: 



2(1 - z)(l - 9z)z^— - z[9z^(d-4) + 10z((i-2)+8-3d] 



dz"^ dz 



+ ((i-3)[z(d + 4) + d- 4] = 12^2-^-6) r2 (^2-0. (2.^ 



The differential equation 1)2. 8|) will be used in Sec. 4 to find the momentum dependence of arbitrary 
periodic constants in the solution of the difference equation (|2.4j) . 
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3 Solution of the dimensional recurrency 

Equation H2.4|) is a second order inhomogeneous equation with polynomial coefficients in d. The full 
solution of this equation is given by (see Ref . and references therein) : 

4''> = 4f + wa{d)4^ + Md)4i\ (3.9) 

where J^p* is a particular solution of H2.4() , Jg^^ , Jg^^ is a fundamental system of solutions of the 
associated homogeneous equation and Wa{d), Wh{d) are arbitrary periodic functions of d satisfying 
relations: 

Wa{d + 2)=wa{d), Wb{d + 2) = Wb{d). (3.10) 

The order of the polynomials in d of the associated homogeneous difference equation can be reduced 
by making the substitution 

j{d) _ ^ (^) -jjd) , . 



r(f -3)r(^ 



The associated homogeneous equation for takes the simpler form 



J 



{d+i) 



27m8(l-z)2(l-9z)2 3 

2d(l-3z)(l-42z + 9z2)z (d+2) (3rf-2)(3d-4) (rf) _ 
27m4(l-z)2(l-9^)2 -^3 3g ^3 - U. 

Putting 

d = 2k-2e, y^''^ = p-'' (3.13) 
we transform Ea. H3.12() to a standard form 

yl^2y(fc+2) ^^B + C k)py'^^+^'^ -{a + k){(3 + k)y^^^ = 0, (3.14) 

where 

16^3 4e(l-3z)(l-42z + 9z2)z 

A = _ -TTT TTT, B - 



27m8(l - z)2(l - 9z)2' 27 m^{l - zf{l - ^z) 
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C = - — , a = -£--, /3 = -£--, (3.15) 

and p is for the time being, an arbitrary constant. In order to get Ea. (|3.14l) into a more convenient 
form, we will define three parameters p, x and 7 by the equations 

Ap"^ =x{l-x), Bp = 'J - {a + P + l)x, Cp = l-2x. (3.16) 

These have the solution 

1 - 2Cp (1 - 9Z)2 ^ g2(^2_ 9^2)2 

2 (l + 3z)3 (g2 + 3m2)3 ' ^'^■^'> 



1 27m4(l-z)2(l-9z)2 27 m2(g2 _ ^2)2(^2 _ 9^2)2 



(3.18) 



V^A + C72 4 z(l + 3z)3 4 (g2 + 3m2)3 

7 = 5p+ + + l)x = -e, (3.19) 
and Eq. 1)3. 14() can accordingly be written in the form 

x(l - x)2/(''+2) + [(1 - 2x)k + 7 - (a + /? + l)x]y(^+^) - (q + A:)(^ + k)y^'''> = 0. (3.20) 
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The fundamental system of solutions of this equation consists of two hypergeometric functions |25j . 
For example, in the case when |1 — x| < 1 (large q^) the solutions are 



(fc) 



(k) 

2/2 



, T{a + k)T{P + k) 
r(a + /3-7 + fc) 



2F1 (a + fc, /3 + fc, a + /3 - 7 + A: + 1; 1 



(1 - x)'' 



2^1(7 - a,7 - /?, 7 - a - /? + 1 - /c; 1 - x). 



(3.21) 



Once we know the solutions of the homogeneous equation, a particular solution Jg^^ can be obtained by 
using Lagrange's method of variation of parameters. Lagrange method for finding a particular solution 
is well described in [JS^. The application of the method is straightforward but tedious. Explicit result 
will be given in the next section. 

It is interesting to note that the argument of the Gauss' hypergeometric function is related to the 
maximum of the Kibble cubic form |26j : 



$(s, t, u) =stu-{s + t + u)rr?{rr? + c^) + Im^irr? + 3g^), (3.22) 
provided that the following condition is satisfied: 

(3.23) 

= I (g^ + 3m^) and we see that the 



s + t + n = g^ + 3m^. 



The maximal value $1, 



9m 



2\2 



occurs at s = t 



kinematical variable 1)3.171 ) can be written as 



X 



stu 



=f=«=i (g2+3m2) ■ 



(3.24) 



This observation may be useful in finding the characteristic variable in the general mass case [Uj . Also 
one can try to apply the method described above to find the imaginary part of the sunrise integral in 
the general mass case in arbitrary space-time dimension. 



4 Explicit analytic expression for 

To find the full solution of Eq. (|2.4() we assume that q"^ is large. The region of large momentum 
squared corresponds to a; ~ 1 and therefore as a fundamental system of solutions of the homogeneous 
equation we take y[''^ and y2^^ . According to (|3.11|1 , (|3.13|1 and (|3.21|1 the solution of the associated 
homogeneous difference equation will be of the form 



J. 



(d) 



3,h 



+ W2{z] 



r(f-i)r(i-|)r(^) 
r(f)r(f-3)r(^) 



d-d 



1 d 

3' 2 
d . 

2 ' 



2 . 

3 ' 1 



p2(d_2 



d 
p2 



r(f -3)r(^bi) (1 



.T2 

X 2 



Fi 



1 2 

3' 3. 



2 ' 



(4.25) 



The arbitrary periodic functions wi{z) and W2{z) can be determined either from the d ^ 00 asymp- 
totics or using the differential equation (|2.8|) . Substituting ()4.25|) into ()2.8j) we obtain two simple 
equations 



z{l - z)(l + 3z)(l - 9z 



^ d'Wi{z) 
dz 



z(\ + 3z)(l - 9z) — + 3(1 - z)w2{z) 
dz 



2(1 + 62-39^2)^1(2) =0, 
0. 



(4.26) 
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Both equations are independent of d and their solutions 

Kiz2(l + 3z)2 



wi{z) 



(l-9z)2(l-z)2 



, W2{z) 



K2Z 



:i + 3z)(l-9z)2' 



(4.27) 



determine the periodic functions up to integration constants ki, K2 which we fix from the first two 
terms of the large momentum expansion of jjf"^ presented in 



r(-i + 2£)r3(i-e)^ 6r2(-e) 



zr2(l + e)r(3-3e) ' ~' ' r(3 - 2e) ' 
The application of Lagrange's method of finding a particular solution gives 



+ 0{z) 



J. 



2d-6 



3p 



(1 + V^Y 



d 



where 



R 



I d-l d 



-,d-l;^/^R,R 



4^^ 



F2 is the Appell function |2H] defined as 

00 

F2(a,/3,/3',7,7';x,?/) = ^ 

k,l=0 



{a)k+0)k{fi')ix^ y' 



k\ i\ 



\x\ + \y\ < 1, 



(4.28) 

(4.29) 
(4.30) 

(4.31) 



and {a)n = T(a + n)/T(a) denotes the Pochhammer symbol. Collecting all contributions, setting 
d = A — 2e, applying Euler's transformation for the first 2F1 function in 1)4. 25() we obtain the following 
solution of the difference equation (|2.4j) : 



J. 



6r2(-e)r2(l+e)(-z)^(l-z) 
r(3-2e)(l + 3z) 



2-2e 



iFl 



27(1 -z)2z 



2-e- (l + 3z)= 



27(1 - zfz 
e; (l + 3z)3 



r(-l + 2£)r3(l - e){-zf'{l - 9zf~^' 
^ m4^-2 r(3 - 3e)z(l + 3z) ^ ^ 

Qr7r)2-4e / 1 3 

+ ^Y^y=^r2 (-1 + e) F2 (^1, 2 - 3 - 2e; ^^i?, 



(4.32) 



This is our main result. 



The imaginary part of J^f^ on the cut comes from the first two terms in 1)4. 32|) . The analytic 



continuation of the two 2F1 functions gives the relatively simple expression: 



ImJ. 



{d) 



-4z TT^VSnm^-'^' 



r(|-e)r(2-£)(l + 3z) 



108^2 



l-e 



2F1 



3531 (1 9;z) 



2-e; (l + 3z) 



(4.33) 



At c? = 4 for the imaginary part this verifies the result of Expanding j!^'^ at q'^ = 9m? we 
reproduced the singular and finite in e parts of several terms of the on-threshold asymptotic expansion 
presented in Ref. [Tlj . 

We found the following integral representation for Appell's F2 function in 1)4.32(1 



F2{l,^,^-e,2-e,3-2e;^R,R 



2T{3 - 2e) 



dt [t(l -t)]2-^ 
{Azt + 1-Z + L) 



2^1 



l,e; Atz + l-z-L 
2-e; 4tz + l- z + L 



, (4.34) 
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where 



L 



4^ 



1 



4z 



(4.35) 



This integral representation can be used for the e expansion of the F2 function. Integral representation 
for the 2F1 functions convenient for e expansion is given in the Appendix. 

The analytic continuation of J^f^ valid near the singular points = 0, m^,9m^ can be directly 
obtained from 1)4. 32() by performing the analytic continuation of the hypergeometric functions. Explicit 
formulae of the analytic continuations of j'^^ in terms of Olsson functions [20] as well as the e expansion 
of the result we are planning to present in a separate publication j3Uj . 

Using (|4.32)) we can find the on-threshold value of the integral. In formula (|4.32j) at = 9m^ 
the imaginary part of the first term cancels the imaginary part of the second term. The real part 
coming from the 2F1 terms cancels the term with 2-^1 which comes from the Appell function F2 at 
the threshold 



Fo 1 



+ 



' 2 

16(1 



1 

e, -,3 



2' 

e) 



2e,2 



3 1 

-'4' 4 



3(1 - 2£) 



3F2 



1,-1 + 2e 
i + e,2 



' 2 



4r(3 - 2e)r(2 - £)T (-i + g) r(2 - 2e) 
H s -^7 ^ 7Z \ 2F1 



3-2-' r2 I 



0r( 



2e 



+ e,2-2e: 



2e: 



(4.36) 



The cancellation of the 2F1 functions happens due to the fact that 



2F1 



2-2£,-i+£; 
- — 2f • 

2 ' 



2l-^-r(f-g)r(|-g) 
3l-r(|-e)r(|-.) 



(4.37) 



Adding contributions from different hypergeometric functions gives a rather simple expression for the 
diagram at q^ = 9m^ 



J. 



id) 



VHe) 



[l-e){l-2e) 



3F2 



1,-1 + 2e, 



1 



2 + e, 2 - e ; 

3 9 75 Svr 
2^2 + 4^ + Y"7! 



+ 0(e)}. 



(l_e)(l-2e) 

The first several terms in the e expansion are in agreement with the result of Ref . ^1] . 



(4.38) 



5 Conclusions 

Here we would like to add several remarks, which underline the most important points which follow 
from the results of this paper. 

Using a new method, for the first time, we were able to obtain an analytic expression for the 
two- loop sunrise, self-energy diagram with equal mass propagators. Until now all attempts to find 
such a result with other methods failed for this integral. This clearly demonstrates that the method 
of dimensional recurrences is a powerful tool for calculating Feynman integrals. In our opinion there 
is a deep reason why, for example, with the differential equation method ^Hj used in an explicit 
formula could not be found. It turns out that even the associated homogeneous differential equation is 
rather complicated. Relevant for this case is the Heun equation ffl,_ with four regular singular points, 
located at = 0,'m'^,9m?,oo. In general the reduction of the Heun equation to the hypergeometric 
equation is a complicated mathematical problem |32j which is not completely solved until now. 
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At the same time the associated homogeneous difference equation for J3 ' is rather simple, and 
admits a reduction to a hyper geometric type of equation with hnear coefficients. 

In fact this is a rather general situation. Kinematical singularities of Feynman integrals are located 
on complicated manifolds. In the case when the differential equations are of the first order there are 
no problems to solve them. However, to solve a second or higher order differential equations in general 
will be a problem because of complicated structure of the kinematical singularities. 

The location of the singularities of Feynman integrals with respect to the space time dimension d 
is well known. This has been used for a rather evident rescaling of the integral by ratios of F functions 
which allowed us to significantly reduce the order of the polynomial coefficients in the difference 
equation as we have seen in Sec. 3. Finally this simplification allowed us to obtain the explicit result. 

We expect that a further development of the method we used in the present paper will help to find 
analytic results for other more complicated Feynman integrals. 



6 Appendix 

The series representation for the Appell function F2: 

F,ia, /?, /?', 7, 7'; X, y) = f; ^-TTTT + ^' ^^'^ ^'^ (6-39) 

is convenient for the analytic continuations and also for the evaluation of F2 at some particular values 
of their arguments. The most frequently used integral representations for F2 is 

F2(a, /?,/?', 7, 7>,^)- r(7)F(7') 



F(/3)F(/3')F(7-/3)F(7'-/?') 



''\udvuf^-^v'^'-\l-uy-'^-\l-vy'-'^'-Hl-uw-vzy. (6.40) 



/O JO 

For the special parameters which appear in the explicit result 1)4. 32() we have: 

/ 3 1 

F2 (1, - - e, -,3 - 2e,2 - 

r(3 - 2e)F(2 - e) dudv [u{l - u){l - v)]^-"^ 



F3 (i - r (i) -'o Jo {l-uw- vz) 



Furthermore, we found the following specific relation between the 2F1 function and yet another 
Appell's function Fi 



2F1 



27(1 -zfz 
2-e; (l + 3z)3 



(l + 3z) ^ /I 1 1 iz 42 , , , , 
^ ^-Fi -, --+£,--+£, 2-e;-— — . (6.42) 



(1-^) ' V2' 2 ' ' 2 ' ' ' (1 + V^)2' (1- V^)2 



To our knowledge such a relation has not been found sofar in the mathematical literature. The integral 
representation for the Fi function reads 

/II 1 \ r(2 - e) du xnl . . ^ 

Fi(-,-- + e,--+e,2-e;w,z)= ^ > [(1 _ n)(l - t«n)(l - zn)] (6.43) 

V2 2 2 y r j F (I - ej -'0 

and is convenient for the e expansion. This -Fi function can be considered as a generating function of 
a new generalization of elliptic integrals which may appear in evaluating Feynman integrals. 
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